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ABSTRACT 

We present a comparison of the properties of substructure halos {subhalos) orbiting within 
host halos that form in Cold Dark Matter (CDM) and Warm Dark Matter (WDM) cosmologies. 
Our study focuses on selected properties of these subhalos, namely their anisotropic spatial 
distribution within the hosts; the existence of a "backsplash" population; the age-distance 
relation; the degree to which they suffer mass loss; and the distribution of relative (infall) 
velocities with respect to the hosts. We find that the number density of subhalos in our WDM 
model is suppressed relative to that in the CDM model, as we would expect. Interestingly, our 
analysis reveals that backsplash subhalos exist in both the WDM and CDM models. Indeed, 
there are no statistically significant differences between the spatial distributions of subhalos 
in the CDM and WDM models. There is evidence that subhalos in the WDM model suffer 
enhanced mass loss relative to their counterparts in the CDM model, reflecting their lower 
central densities. We note also a tendency for the (infall) velocities of subhalos in the WDM 
model to be higher than in the CDM model. Nevertheless, we conclude that observational tests 
based on either the spatial distribution or the kinematics of the subhalo population are unlikely 
to help us to differentiate between the CDM model and our adopted WDM model. 

Key words: methods: n-body simulations - galaxies: halos - galaxies: evolution - cosmology: 
theory - dark matter 



1 INTRODUCTION 

The currently favoured ACDM model of cosmological structure 
formation has proven to be extremely successful at describing the 
clustering of matter on intermediate to large scales (e.g., Springel 
|et al.|2005l|Spergel|2007l >. In contrast, it has been argued that the 
predictions of the ACDM model are at odds with observations on 
the scales of galaxies, on the basis of cosmological iV-body simula- 
tions. Cold Dark Matter (CDM) halos are predicted to have "cuspy" 
density profiles with inner logarithmic slopes of approximately -1.2 
(e.g., |Navarro et al.|20 04, Tasitsiomi et al. 2004; Diemand et al. 
|2005||Reed et al.|2005| l, whereas high resolution observations of 
low surface brightness galaxies appear to require halos with con- 
stant density cores (e.g., |Gentile et al.|2007||McGaugh et al.|2007) . 
Furthermore, CDM halos are predicted to contain a wealth of sub- 
structure, which we might expect to observe as satellite galaxies 
within galactic halos, in sharp contrast to the observed abundance of 
satellites around our Galaxy and others ( Kl ypin et al.| 1 999 , Moore 
|etal.|1999| >. 

Suggested solutions to these problems have included allowing 
the dark matter to be collisional (i.e. self-interacting) rather than 
collisionless < |Spergel & Steinhardt| [2500} |Bento et aL]|2"0"00"l l, 



allowing it to be warm rather than c old ([Bode et alT||2001| 
|Avila-Reese et al.|[200Tj |Knebe et al.||2002) , and introducing 
non-standard modifications to an otherwise unperturbed CDM 
power spectrum (e.g., Bullock 2001 , Little et al. 2003). Arguably 
the most promising (and least intrusive) modification to the dark 
matter paradigm is to allow the dark matter particle to be warm. 
In such a case, warm dark matter particles will have a relatively 
high thermal velocity dispersion at decoupling and therefore a 
non-negligible free-streaming scale A&. This modification results in 
a change to the primordial matter power spectrum, corresponding to 
a damping of density perturbations on scales below & filtering scale 
Rf (which is related to the free-streaming scale Af s ), which in turn 



is related to a filtering mass Mf i Bardeen et al. 


|1986||Bodeetal.| 


|2001 ||Avila-Reese et al.|2001||Knebe et al.|2002 





Previous studies have revealed that WDM can resolve some 
of the tension between theoretical prediction and observation. In 
particular, the abundance of substructure halos (hereafter subhalos) 
is greatly reduced in WDM models (Bode et al. 2001 , Avila-Reese 
|et al.|200l[|Knebe et al.|2002^ compared to the CDM model. How- 
ever, the simulations used in these studies did not have sufficient 
resolution to follow the orbits of subhalos in detail, and so these 
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results are based on a snapshot of the subhalo population, static in 
time rather than a dynamic entity. 

Over the last decade cosmological iV-body simulations have 
advanced and reached a stage where it is possible to study the dy- 
namics of well resolved subhalos and satellite galaxies and use them 
as a probe of cosmology. This field - dubbed "near-field cosmology" 
(e.g., Bland-Hawthorn & Peebles 2006 ) - has prompted numerous 



studies (e.g. 



Warnick et al. 2008, Colin et al. 2007, Sales et al 



20071 [Faltenbac her et al.|2007||Diemand et al.||2007| |Libeskind| 



t al.| 
rndj 



et al.|20 07 , Kang et al. 2007 ; War nick & Kne be 2006 , Agust sson &| 
Brair.erd|20061|Knebe et al.|2006||Libeskind et al.|2005||Reed et"ak| 
20051 |De Lucia et al.|2004l|Gao et al.|2004l|Gill et al.|20041|Knebel 
et al.|2004[|Kravtsov et aTf2 004; Zent ner & Bullock|2003) , but as 
yet there have been no in-depth comparisons of subhalos in CDM 
and WDM models. We address this in the present study, using high 
resolution resimulations of a set of "identical" galaxy clusters (see 
next section) forming in CDM and WDM models to compare and 
contrast properties of their subhalo populations. 

Rather than seeking to reproduce and verify all of the recent 
results for subhalos derived from simulations of the CDM model, 
we focus on selected properties of the subhalo population. Namely, 
we set out to validate the existence of the so-called "backsplash" 



population reported in 


Gill et al. ( 2005} (see also |Mamon et al. 


2004||Moore et al.|2004 


Balogh et al. 2000), the anisotropic spatial 


distribution of satellites (e.g., Bailin et al. 2007, Agustsson &| 


Brainerd 2006 , Zentner 


et al. 2005 , Libeskind et al. 2005 , Knebe 



et al. 20041, the putative age-distance relation (cf., |Gao et al.|2004| 
Moore et al. 2004 1, and the degree to which subhalos suffer mass 
loss. Wc also examine the relative velocities of subhalos with respect 
to their hosts. This is of particular interest because it has been argued 
that the collision velocity of the "Bullet Cluster" (Milosavljevic 
|et al.|2007||Markevitch et al.|2002|> may pose another challenge 
to CDM (e.g., |Springel & Farrar|2007l |Hayashi & White|2006l >. 
While it might be expected that the large scale tidal field will be 
more important for collision velocities in mergers, it is nevertheless 
interesting to see whether the precise nature of the dark matter 
might play a role. Having characterised these properties, we can 
quantitatively address the question of whether or not the spatial and 
kinematic properties of subhalos at the present day could be used to 
differentiate between the CDM and WDM models. 



2 THE NUMERICAL SIMULATIONS 
2.1 The Raw Data 

Our analysis is based on suite of high-resolution iV-body simula- 
tions. They were carried out using the publicly available adaptive 
mesh refinement code MLAPM ( |Knebe et al.|200l) focusing on the 
formation and evolution of dark matter galaxy clusters containing of 
order one million particles, with mass resolution 1.6 x 10 8 h' 1 M Q 
and force resolution ~2/i _1 kpc. We first created a set of four inde- 
pendent initial conditions at redshift z — 45 in a standard ACDM 
cosmology (Q = 0.3, fix = 0.7, fi 6 = 0.04, h = 0.7, cr 8 = 0.9). 
512 3 particles were placed in a box of side length 64/i _1 Mpc giv- 
ing a mass resolution of m p = 1.6 x 10 8 ft _1 Mq. For each of 
these initial conditions we iteratively collapsed eight adjacent parti- 
cles to one particle reducing our particle number to 128 3 particles. 
These lower mass resolution initial conditions were then evolved 
until z — 0. Then, as described by Tormen ( 1997 1, for each cluster 
the particles within five times the virial radius were tracked back 
to their Lagrangian positions at the initial redshift (z = 45). Those 



particles were then regenerated to their original mass resolution and 
positions, with the next layer of surrounding large particles regen- 
erated only to one level (i.e. 8 times the original mass resolution), 
and the remaining particles were left 64 times more massive than 
the particles resident with the host cluster. This conservative crite- 
rion was selected in order to minimise contamination of the final 
high-resolution halos with massive particles. 

The three warm dark matter halos were simulated using the 
same techniques. In fact, the only difference between the CDM 
and WDM halos is the functional form of the primordial power 
spectrum used as an input for the initial conditions generator. We 
follow Bar deen et al.| fl986l and modify the CDM power spectrum 
by multiplying it with a damping function F^ f - DM (k), where 



-FwDM(fc) = exp 



kR d 
2 



{kR d y 
2 



(1) 



Following Barde enet al.|{1986") , we parameterise the damping scale 
Rd in terms of the warmon density parameter f2 w dm, its mass m w dm 
and the dimensionless Hubble parameter h, 



R d ^ 0.074 (^) 



-4/3 fa 



0.3 



1/3 



07 



5/3 



/i _1 Mpc.(2) 



We adopt a warmon mass of m w( j m = 0.5 keV, which gives a 
damping scale of R d =0.186/i~ 1 Mpc. The filtering scale Rf 
can be obtained by determining the wavenumber of the mode at 
which the amplitude of the linear density fluctuation is suppressed 
by a factor of two, and then computing half the comoving wave- 
length. It is straightforward to evaluate this from the WDM and 
CDM power spectra. For our choice of warmon mass and cosmo- 
logical parameters, we find that this wavenumber corresponds to 
k = 2.553/i Mpc -1 and therefore the filtering scale (mass) is 
R f = 1.24/i" 1 Mpc(M / = 6.653 x lO 11 /i _1 M ). 

Note that our choice of warmon mass is lower than recent lower 
limits derived from combined analysis of observed properties of the 
matter power spectrum as inferred from the Sloan Digital Sky Survey 
Lyman-a flux power spectrum, cosmic microwave background data 
and the galaxy power spectrum, which vary between, e.g., m w d m J? 
3 keV l |Abazajian|2006"} and m wdm > 10 keV ( |Viel et al. 2006). 
However, it is consistent with published estimates of the ~ 0.5 keV 
warmon mass that would resolve the "overabundance" of dark matter 
substructure in galactic halos (e.g. |Moore et al.| F9 99 Dalca nton &| 
Hogan 2001 , Goerdt et al. 2006). By focusing on the lowest warmon 
mass that could be considered consistent with observational data, 
we can explore subhalo dynamics in a plausible model, yet one in 
which the effects of the warmon should be more pronounced and 
therefore easier to identify when comparing and contrasting with 
the CDM model. 

We note that the three CDM halos CDM1, CDM2, and CDM3 
have appeared previously in both Warnick & Knebe ( 2006} and 
| Warnick et a l. (2008), in which they corresponded to the "C3", "C7" 
and "C8" systems. 



2.2 Discreteness Effects 



It has been argued in the recent study of Wang & White ( 2007} that 
WDM halos below a given fraction of the filtering mass Mf are 
spurious, arising from the unphysical fragmentation of filaments. 
They provided the following expression (based upon simulations of 
the Hot Dark Matter model) 



Mum ~ 10.1pdfc pe ; 



(3) 
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~p is the mean density, d is the mean interparticle separation, and 
fcpcak is the wavenumber at which A 2 (fc) = k 3 P(k) reaches its 
maximum. 

We have taken care to compute Mn m , noting that our simula- 
tions use boxes of side 64/i _1 Mpc, an effective number of particles 
512 3 , a density parameter of f2o=0.3 and a peak wavenumber of 
fc pcak =1.78 h Mpc _1 Q These numbers give d=0.125 ft _1 Mpc, 
p=8.3265 h 2 M Q Mpc -3 and so M iim =3. 317828 x 10 10 /i _1 M . 
Our particle mass is m p = 1.626269 x 10 8 /i _1 M Q , and so we find 
that Mii m is equivalent to 204 particles. Therefore, this corresponds 
to the mass cut applied in the following analysis. 

We have also applied a cut of 2M]i m to our data and checked 
our results to ensure that they are unaffected by spurious halos. This 
reveals that our results remain unchanged and therefore are stable. 
We do not find any systematics biases in our data if we employ mass 
cuts of Mii m or 2Mii m . We have looked for trends in our data that 
we would expect to be present if they were affected by spurious 
haloes (such as the distribution of concentrations) but we find no 
obvious signatures. Therefore we conclude that our results are robust 
and unaffected by particle discreteness. 



2.3 The Halos 

Both the halos and their subhalos are identi fied using AHEpj a modi 



Gill et al 



(20041, which 



fication of the MHEQalgorithm presented in I 
has been parallelised using the MPI (Message Passing Interface) li- 
braries. AHF utilises the adaptive grid hierarchy of MLAPM to locate 
(sub)halos as peaks in an adaptively smoothed density field. Local 
potential minima are computed for each peak and the set of particles 
that are gravitationally bound to the peak are returned. If the peak 
contains in excess of 20 particles, then it is considered a (sub)halo 
and it is retained for further analysis. 

For each (sub)halo we calculate a suite of canonical properties 
from particles within the virial/truncation radius. We define the virial 
radius i? v i r as the point at which the density profile (measured in 
terms of the cosmological background density pb) drops below the 
virial overdensity A vir , i.e. M(< i? vir )/(47ri? 3 ir /3) = A vir p b - 
Here pb is the mean density of the background (Universe). Follow- 
ing convention, we assume the cosmology- and redshift-dependent 
definition of A v i r ; for a distinct (i.e. host) halo in a ACDM cos- 
mology with the cosmological parameters that we have adopted, 
A v i r = 340 at z — 0. This prescription is not appropriate for subha- 
los in the dense environs of their host halo, where the local density 
exceeds A v i r pb, and so the density profile will show a characteristic 
upturn at a radius R <J R V1I . In this case we use the radius at which 
the density profile shows this upturn to define the truncation radius 
for the subhalo. Further details of this approach (and especially the 
"halo tracking" used to obtain the temporal evolution of subhalos) 
i |Gill et al.|p004l >. 



can be found in 



In Table [T we summarise some of the properties of the host 
halos along with particulars of their respective subhalo populations. 



The peak value fc pe ak = l-78 h Mpc -1 has been determined numerically 
from the (tabulated) warm dark matter power spectrum used in this study 
based upon a CDM power spectrum calculated with the publically available 
CMBFAST code (Seljak & Zaldarriaga 19961 and modified according to 

AMIGA's-Halo-Finder; AHF can be downloaded from 

http://www.aip.de/People/aknebe/AMIGA. AMIGA is the 
successor to MLAPM. 
2 MLAPM's-Halo-Finder 
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Figure 1. Minimum distance as a function of present-day distance. 



3 ANALYSIS OF THE SUBHALO POPULATION 
3.1 The "Backsplash" Population 

It has been noted that a significant population of halos on the 
outskirts of present-day galaxy- and cluster-mass host halos once 
resided within the virial radii of these hosts at earlier times jWarnick| 



etal.|20 08; Gill e t al.|2005||Mamon et al.|2004||Moore et al.|2004| 
Balogh et al. 2000 ). These results are based on simulations of the 
CDM model, but it is interesting to ask whether such a "backsplash" 
population exists in the WDM model. We expect there to be fewer 
satellites in WDM models and these satellites will tend to have lower 
concentrations (Knebe et al. 2002; Avila-Reese et al. 2001; [Bodel 
|et al.|20 01 1, which, when combined, should affect the numbers of 
"backsplash" halos. The lower the concentration of a (sub)halo, the 
greater the likelihood that it will be tidally disrupted within the host. 
Therefore, we might expect the numbers of backsplash halos to be 
suppressed in WDM models. 

In Figure[T]we plot the minimum halocentric distance reached 
by a (sub)halo against its present day halocentric radius. Note that 
we combine data for all three halos in the CDM and WDM models 
respectively. As expected, backsplash halos are present in both the 
CDM and WDM models, although the numbers are reduced in the 
WDM model. Table[T]reveals that the youngest system (i.e. host #3) 
has the smallest fraction of backsplash halos; this reflects the fact 
that this system has experienced a recent triple merger (cf. Warnick 
|et al.|2008[|Warnick & Knebe|2006[ l. Interestingly, we find again 
(cf. |Gill et al.|2005| > that the number of infalling halos is of the 
same order as the number of backsplash halos in both the CDM and 
WDM models. 
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Table 1. Summary of the host halos properties and their subhalo populations. The age is given in Gyrs, i? v ir is measured in /i" 1 kpc, masses in 10 /i Mq 
and the velocity dispersion a v in km/sec. We follow Lacey & Cole ( 1 993 1 and define the formation redshift of our host halos as the redshift at which the 
halo's most massive progenitor first contains in excess of half its present day mass. The concentration Ci /g = R vlr / Ri /g is defined via the radius R 1 /g that 
encompasses l/5 th of the virial mass. Shape is quantified by the triaxiality parameter T = (a 2 — b 2 )/ (a 2 — c 2 ) and the eigenvalues of the inertia tensor 
a > b > c. A^ t measures the number of a certain subset X (int=interior, inf=inf ailing, back=backsplash) of subhalos while N^?' 5 " gives the total 
number of subhalos within 2.5i? v j r ; note that these numbers only reflect subhaloes in excess of 200 particles. 
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Figure 2. The cumulative distribution of all interior subhalos normalised 
to the total number of satellites. The short thin lines represent curves with 
the logarithmic slopes of a = —0.8 (solid) and a = —0.6 (dashed), 
respectively. 



3.2 Mass Spectra 

The mass spectrum of satellite galaxies in WDM and CDM has 
been studied previously l |Rnebe et al.|2002||Avila-Reese et al.|2001| 
|Bode et alj2 001 ), but we consider it here briefly for completeness. 
In Figure [2] we show the cumulative mass functions of subhalos 
normalised to the total number of subhalos, where we normalise the 
subhalo mass by the host halo mass. It is readily apparent that the 
abundance of low-mass halos is suppressed in WDM cosmologies. 
We fit a power-law to this mass function, 



N(> M) 



total, interior 



Ms, 



Mvir,l 



(4) 



and obtain logarithmic slopes of a — —0.8 for CDM and a = —0.6 
for WDM. This is consistent with the findings of previous studies 



for CDM subhalos (e.g., Shaw et al.|2007 


Reed et al. 2005||Gill| 


|et al.|2004||De Lucia et al. 2004[|Gao et al. 


2004|. 



We also calculate the mass functions of the infalling and back- 
splash halo populations, which we show in Figure[3] We observe a 
general trend that backsplash halos contain fewer high-mass objects 
in comparison to the infalling halos. This reflects the importance of 
tidally induced mass loss for backsplash halos, which we quantify 
in the next section (cf. also Wamick et al. 2008). 
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Figure 3. As in Figure [2] but for the infalling and backsplash populations 
respectively. 



3.3 Mass Loss 

There is a general consensus that subhalos in WDM models are less 
concentrated than their counterparts in CDM models ( |Knebe et al.| 
|2002[|Bode et al.|2001[|Avila-Reese et al.|200l| and therefore more 
susceptible to tidal destruction while orbiting within the dense envi- 
rons of their host halo. We verify this in Figure[4]and Figure|5]where 
we plot the total (fractional) mass loss as a function of distance to 
the host for both the interior and the backsplash population. The 
mass loss is measured over the time period from infall onto the host 
(i.e. the first time a satellite crosses the virial radius of the host on 
an inward trajectory) until the present-day. 



Subhalos in WDM Models 5 



1.0 



a 0.8 



1 0.6 



0.4 



0.2 



3x CDM 



0.0 
1.0 



2 0.8 



1 0.6 



0.4 



0.0 



# ++ + + + + " * + + 
+ > + + : + + ; ' ? i J 

: I H I t Y + 1 4v 



3x WDM 



0.0 0.2 0.4 0.6 O.i 

D mm/ R »ir,tiost 



1.0 



S 0.8 



1 0.6 



0.4 



0.2 



0.0 
i .0 



a 0.8 



1 0.6 



0.1 



0.2 



0.0 



+ 

+ + + + 

+ 


4- 

++ 


' , + + 
+ ^^ 

1 — 


+ + + + 

+ ^ 

4- 

+ 


4- 

4- 



3x CDM 



— * 



3x WDM 



0.0 0.2 0.4 0.6 O.i 

Dm ln /R rt r.ho.t 



1.0 



Figure 4. The total (fractional) mass loss of interior subhalos as a function 
of minimum distance to the host. The error bars represent the standard 
deviation. 

In both dark matter models the average mass loss (presented 
as histograms in Figure [4} is a monotonic decreasing function of 
minimum distance. However, in the WDM model this function is 
pointwise greater than the corresponding curve in CDM, which 
would be expected if mass loss is enhanced as a result of the lower 
concentrations of subhalos in the WDM model. 

Surprisingly, the relation between mass loss and minimum dis- 
tance is not as steep for backsplash halos as for interior subhalos 
(Figure [5J. Nevertheless, (sub)halos plunging deeper into the po- 
tential well of the host experience greater mass loss - as expected 
and confirmed for the interior population ( e.g.,|Sales et al.|2007| 
IDiemand et al.|2007||Kravtsov et al.|2004||Gao"et al.|2004| l. On 
average, backsplash halos in the WDM model suffer greater mass 
loss than their CDM counterparts. 

3.4 Spatial Anisotropy 

There is good reason to believe that the spatial distribution of sub- 
halos (and the subset corresponding to satellite galaxies) in both 
cluster sized systems and galactic halos is anisotropic in the CDM 
model (e.g., Faltenbacher et al. 2007 ; Kang et al. 2007) |Libeskind| 
et al.|2007||Agustsson & Brainer d|2007| |2006| |Zentner et al. |2005 1 
Libeskind et al.|2005[|Knebe et al.|2004| >, and so it is interesting to 
ask whether the same can be said of the WDM model. 

For each of our host halos, we compute the cumulative fraction 
of subhalos which have cosine of the angle 

cos$ = R sat ,z=o ■ Ei; (5) 



Figure 5. Same as Figure[4]but for the backsplash population. 



this measures the position of a subhalo relative to the host R sa t,z=o 
and the host's major axis Ei. The host's major axis is identified 
using the eigenvalues and eigenvectors of its moment of inertia ten- 
sor, where the eigenvector corresponding to the smallest eigenvalue 
defines the major axis. 

The resulting distributions of cos $ are shown in Figure|6]for 
subhalos within the host's virial radius, which confirms that the 
spatial anisotropy is present in the WDM model although not as 
pronounced as for the CDM case. Although we do not show the 
result here as it (probably) lies below the credibility level of the 
WDM simulation (Wang & White 2007), we note that very low 
mass systems (i.e. M sa t < 10~ Mh os t) correlate more strongly 
with the major axis than the remainder of the satellites - for both 
dark matter models. This is consistent with the expectation that 
(especially low mass) objects are primarily channelled along the 
filaments feeding the cluster (cf. also |Knebe et al.||2004| l. The 
thin dashed line represents an isotropic distribution or the "uniform 
continuous distribution function" (UCDF). 

We have computed the same distribution for both the back- 
splash and infalling satellites in Figure [7] Interestingly we find that 
the spatial anisotropy is even stronger for the backsplash popula- 
tion than for the interior objects even though it is skewed towards 
cos $ ~ 0.55 for the WDM model and hence no perfect alignment 
with Ei anymore (but nevertheless an anisotropic distribution). One 
possible explanation for this could be that backsplash halos tend 
to be on radial orbits that are either plunging through the host or 
grazing the "virial surface". If this is the case, we might expect to 
observe a tendency for the infalling population to align with the 
major axis as verified in the lower panel of Figure [7] In addition 
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Figure 6. Cumulative fraction of satellites with the absolute value of the 
cosine of the zenith angle < | cos</>|. The zenith angle, < <j> < n, is 
defined as the angle from the major axis of the dark matter distribution of the 
host. The dotted line corresponds to an isotropic distribution or the "uniform 
continuous distribution function". 
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Figure 7. Same as Figure|6]but for the backsplash (upper panel) and infalling 
(lower panel) population, respectively. 
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Figure 8. Age-Distance Relation. We define the "age" of a subhalo as 
the period of time that has elapsed since it first crossed the virial radius 
of its host on an inward trajectory. The redshift at which this occurs is 
the "accretion redshift". D z= o/i? v ir host measures the subhalo's current 
distance to with respects to the centre of the host. For clarity, the triangles 
represent backsplash halos while black crosses denote "interior" subhalos. 



we note that the infalling WDM satellites show a more pronounced 
tendency to be aligned with the major axis of the respective host. 
This goes along with the expectation for (sub)haloes to be concen- 
trated in the filaments and channelled along them into the cluster 
(e.g., |Knebe et al.|2004) . 



3.5 Age-Distance Relation 

Is there is a correlation between the infall time of a satellite and its 
present-day halocentric radius? Previous studies using cosmological 
simulations have argued that "older" subhalos tend to lie closer to 
the centre of the host ( |Willman et al.|2004[ |Gao et al.|2004) , but 
there are also claims to the contrary (Moore et al. 2004). Here the 
age of a subhalo corresponds to the period of time that has elapsed 
since it was first accreted by the host (accretion redshift), entering 
the virial radius on an inward trajectory. 

We investigate whether such a correlation between a subhalo's 
age or accretion redshift and halocentric radius exists in our data 
in Figure[S] Crosses represent "interior" subhalos within the virial 
radius at the present day while triangles represent backsplash halos. 
Because we output our snapshots at discrete intervals, the times at 
which subhalos are accreted appear discrete. It is possible to correct 
for this discreteness by interpolating the growth of the virial radius 
and the positions of subhalos between snapshots, but we output 
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snapshots sufficiently frequently to make the uncertainty introduced 
by discreteness negligible. 

Figure [8]reveals that the correlation between radius and age is 
not a straightforward one. The CDM1/WDM1 system is hard to in- 
terpret because there is no strong trend for subhalos within the virial 
radius. However, there are interesting trends in the CDM2/WDM2 
and CDM3/WDM3 that suggest that there may be distinct popula- 
tions following distinct age-distance relations. The most recently 
accreted subhalos, with accretion redshifts z <. 0.1, show the ex- 
pected trend for accretion redshift to increase with decreasing red- 
shift. However, subhalos accreted at z J> 0.1 appear to follow an 
inverse relation, tending to have higher accretion redshifts for larger 
halocentric radii, and this trend continues beyond the virial radius 
into the backsplash population. Finally, we note that the "oldest" 
subhalos do not appear to follow any trend, instead forming a hard 
upper edge in each panel. However, this edge is an artifact of our 
method for tracking subhalos (cf . |Gill et alT[ 2004 1 and corresponds 
to the formation redshift, Zf orm , of the respective host halo. This 
explains why there is a systematic shift to higher redshifts from 
the lower plot to the upper one - our halo tracker starts following 
subhalos at Zf orm and so it cannot "see" (sub)halos prior to Zf orm . 
Therefore all subhalos that resided within the host at this initial time 
appear as infalling ones. 

This figure also reveals that no backsplash halos have been 
accreted more recently than z « 0.15, which corresponds to a 
period of ~ 2 billion years. This might be considered the minimum 
time a backsplash halo spends within the virial radius of the host 
halo. We can compare this to the time scale idyn for a subhalo to 
complete one circular orbit at the virial radius; 




Figure 9. Cumulative distribution of relative velocity between all interior 
subhalos and their respective host. The thin vertical line is representative of 
the collisional speed of the "Bullet" cluster. 



a large scatter. We note that these authors identify subhalos at 2=0 
and track their merger trees backwards in time. 

It is clear from our analysis that any conclusions we draw must 
be tentative - if we are to gain greater insight into the age-distance 
relation then we must draw upon a larger (i.e. statistical) sample of 
host halos. Nevertheless we note that any age-distance relation is 
apparent in the subhalo populations in both the CDM and WDM 
models. 



^dyn — 



where we used 



AvirPb- This leads to t^ y 



(6) 



6 x 10 9 yr. 



Therefore, the minimum time ti ns id c a subhalo spends inside its 
host is approximately 1 /3 rd the time it would take to complete one 
orbit at the virial radius. This suggests that backsplash haloes are on 
preferentially radial orbits, and explains why subhalos accreted at 
earlier times are preferably found outside 7? v i r today. 

We conclude that an age-distance relation akin to the one re- 
ported by |Gao et al.| (2004) is valid only for "recently" accreted 
subhalos. According to Figure [8] there is a clear correlation between 
a subhalo's age and its distance apparent in CDM2/WDM2 and 
CDM3/WDM3, the two youngest sets of hosts in our sample, and 
this correlation is apparent for objects accreted after z = 0.2. How- 
ever, as we note above, there is some evidence that there may be 
distinct subhalo populations, separated according to their accretion 
redshift, that follow distinct relations and inverse relations. 

It is worth noting that |Gao et aT| {2004 ) observed an age- 
distance relation over a much greater time span ranging from 
z ~ 0.9 down to z ~ 0.3 — 0.4. These authors identify all subhalos 
at z=\ and track them forward in time, which is equivalent to our 
approach (we start tracking subhalos at the formation redshift of the 
host, which varies between z ~ 0.8 for CDM 1 /WDM 1 to z ~ 0.2 
for CDM3/WDM3). We note that accretion in their data appears to 
be complete by z ~ 0.3 (see upper panel of their Figure 15). How- 
ever, it remains unclear why there is no further accretion apparent 
in that plot for smaller redshifts as the fraction of accreted subhalos 
increases at least down to z ~ 0.1 (cf. upper left panel in their 
Figure 12). Moore et al. (2004) deduced from their analysis that any 
age-distance relation present in their data had to be very weak with 



It is notable that there are no counterparts in the WDM1 system 
to the "old" backsplash population we observe in the CDM1 model. 
To better understand why this might so, we checked the distribution 
of infall velocities for subhalos that were accreted at early times in 
each of host (i.e. z > 0.75 for CDM1 and z > 0.80 for WDM1). 
We found that the typical subhalo velocity in the CDM1 run was 
approximately 30% larger than in the WDM1 run. This explains 
why we do not find an "old" backsplash population in the WDM1 
run - the typical infall velocity of a subhalo is too low to allow it 
to escape the host and become a backsplash halo. It is interesting 
to note that this behaviour is peculiar to the CDM1/WDM1 set of 
hosts, but we consider it statistically insignificant. 



3.6 Relative Velocities 

So far we have compared and contrasted the spatial distribution 
and subhalos in the CDM and WDM models, and the degree to 
which these subhalos suffer mass loss. It is also of interest to ask 
whether the kinematics of subhalos in the respective models dif- 
fers. We already learnt that the early infalling population in CDM1 
is marginally faster than their counterpart in WDM1. Hence, if 
there are systematic differences between the models, what might 
the implication be for a system such as the "Bullet Cluster"? This 
is an extremely high velocity merger between two galaxy clus- 
ters ( |Milosavljevic et al.|2007[|Markevitch et al.|2002^ and it has 
prompted discussion as to whether such high relative velocities (of 
order ~ 4500 km/sec) can be accommodated within the CDM model 
(|Nusser|2007[|Angus & McGaugh|2007||Springel & Farrar|2007| 
|Hayashi & White|2006| |Gill et al.|2005| >. We might expect such 
high relative velocities to be sensitive to the large scale gravitational 
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Figure 10. Same as Figure[9]but for the backsplash and infalling population, 
respectively. 

fieldJ^Furthermore, if there are differences in the relative velocity 
distributions in the WDM and CDM models, this could also allow 
limits to be placed on the nature of the dark matter. 

In Figure[9]we plot the cumulative distribution of relative veloc- 
ities (Vs a t — Vhoat) for all interior subhalos, where V sa ,t and Vkost 
are the centre-of-mass velocities of all particles inside the virial 
radius of the subhalo and the host, respectively. Relative velocities 
have been normalised to the circular velocity V-ni of the host at the 
virial radius. If we compute this quantity for the "Bullet Cluster" 
using the estimate of the mass deduced from weak lensing ( |Clowe| 
|et al.|2004 ), the normalised collision speed (V sa t — Vtost)/Krir is 
approximately 1.9 (and shown as a thin vertical line). This figure 
reveals that ~6% of subhalos in the CDM model have normalised 
relative velocities in excess of 1.9, compared to ~10% in the WDM 
model. In other words, the probability of a high-speed encounter is 
greater in the WDM model than in the CDM model. 

We have computed the same distributions for infalling and 
backsplash (sub)halos and the results are shown in Figure[l0] There 
we observe that the infalling population in the WDM model is 
marginally faster than its CDM counterpart. However, we also note 
that the fastest infalling satellitecan always be found in the CDM 

3 We caution the reader that the "Bullet" cluster is a system where the host 
(and the "Bullet" itself) is an order of magnitude more massive than the hosts 
(and satellites) presented in this study. While the existence of the "Bullet" 
cluster may serve as a motivation for the study of relative velocities any 
conclusions drawn from our results are to be extrapolated to the "Bullet" 
system with care. This is especially so because the differences between CDM 
and WDM are less prominent on scales corresponding to the "Bullet" cluster. 



Figure 11. The absolute value of the relative velocities between the 
(sub)halos and their respective host versus their distance at the present time 
is shown. Diamonds represent the infalling population and black crosses 
denote the interior and backsplash populations. The thick solid lines are the 
mean values with error bars showing the standard deviation in each bin. For 
clarity, these are offset for the infalling halos. The squares denote the 20 
most massive subhalos. 



model. Or in other words, in WDM there are more infalling subhalos 
with relative velocities up to about 1.2 x V v i T , but in CDM there 
exists the odd satellite with a velocity as high as ~ 1.5 x V v i T . One 
potential explanation of this may be that subhalos suffer "dynamical 
friction within the filaments". As shown by Knebe et al. ( 2003 1, more 
mass in CDM filaments is found in gravitationally bound objects 
whereas the mass in a WDM filament is more uniformly distributed, 
which may lead to enhanced dynamical friction. Therefore, subhalos 
falling along filaments may have their infall velocities reduced and 
hence we a) do not find exceedingly fast subhalos and b) observe 
an increase in the number of objects in the range 0.75-1.2 V v i T . 
The situation though is different for the backsplash population that 
appears to be slower than its CDM counterpart. 

In Figure [TT] we plot the (unnormalised) relative velocities 
of subhalos as a function of their present-day halocentric radii. 
Interior and backsplash populations are represented by crosses while 
infalling subhalos are represented by diamonds; we also highlight 
the 20 most massive subhalos by green squares. 

There are a number of points to note in this figure. The first 
is that the subhalos with the highest relative velocities are concen- 
trated towards the centre of the host, nestled in its potential well. 
The second is that the most massive subhalos (green squares) are 
not responsible for the high-velocity tail that we observe in Figure|9] 
The third point is that the WDM backsplash population a) does not 
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extend spatially as far out as for its CDM counterpart and b) has 
lower velocities leading to the observed steeper decline of the ve- 
locities with increasing clustercentric distance. The fourth and final 
point is that infalling subhalos are a distinct population kinemati- 
cally, tending to have higher velocities than backsplash galaxies (see 
also |Gill et al.|2003] >. We conclude that the kinematics of subhalos 
is unlikely to allow us to differentiate between the WDM and CDM 
models. 



4 DISCUSSION & CONCLUSIONS 

We have compared and contrasted the properties of subhalos orbiting 
in a set of simulated galaxy cluster hosts in the CDM and WDM 
models. The mass and force resolution of our simulations were 
sufficient to our host halos with ~ 10 6 particles within the virial 
radius at 2=0, and we could follow the orbital evolution of hundreds 
of subhalos in detail using outputs finely spaced in time (At ss 
170 Myrs) from the formation time of the host to the present day. 

Our study has revealed that many of the properties of subhalos 
in the CDM model and the WDM model we have studied are similar. 
Subhalos in both the CDM and WDM models are distributed 
anisotropically with respect to the major axis of their host, and the 
phenomenon of "backsplash" halos is common to both models. 
Other studies have shown that low-mass halos in WDM models 
tend to be less centrally concentrated than their counterparts in the 
CDM model (Colin et al. 2007 1, and this leads to enhanced mass 
loss via tidal stripping for subhalos in WDM models. We find no 
evidence for a well pronounced correlation between the age of 
a subhalo and its present day halocentric radius in either model. 
Interestingly, we find that subhalos in the WDM model are likely to 
have higher (infall) velocities than in the CDM model. 

Our results nevertheless suggest that it is unlikely that the 
spatial distribution and kinematics of subhalos can be used to dif- 
ferentiate between the CDM and WDM models at z=0. It might be 
possible to detect differences at higher redshifts, when the effect 
of the filtering mass is more pronounced (e.g. |Power et al. 2008 ). 
Furthermore, it might be possible to detect differences in the stellar 
populations and star formation histories of the satellite galaxies that 
are hosted by subhalos (e.g. jGao & T heuns 20071, which will be 
sensitive to the mass assembly histories of the (sub)halos. However, 
such measures depend explicitly on the veracity of galaxy formation 
modelling, and so it seems more likely that estimates of the small 
scale power spectrum deduced from the Lyman-a forest (e.g. |Viei| 
|et al.|2007| > may provide stronger constraints. Nevertheless, it is 
important to consider the various strands of observational evidence 
when piecing together the dark matter puzzle. 
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